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ABSTRACT 

Gamma-ray bursts (GRBs) are divided into two classes according to their 
durations. We investigate if the softness of bursts plays a role in the conven- 
tional classification of the objects. We employ the BATSE (Burst and Transient 
Source Experiment) catalog and analyze the duration distributions of different 
groups of GRBs associated with distinct softness. Our analysis reveals that the 
conventional classification of GRBs with the duration of bursts is influenced by 
the softness of the objects. There exits a bimodality in the duration distribution 
of GRBs for each group of bursts and the time position of the dip in the bimodal- 
ity histogram shifts with the softness parameter. Our findings suggest that the 
conventional classification scheme should be modified by separating the two well- 
known populations in different softness groups, which would be more reasonable 
than doing so with a single sample. According to the relation between the dip po- 
sition and the softness parameter, we get an empirical function that can roughly 
set apart the short-hard and long-soft bursts: SP = (0.100 ± 0.028)T 90 ( ' ' 85±0 ' 18 \ 
where SP is the softness parameter adopted in this paper. 
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1. Introduction 

The bimodality in the duration distribution of GRBs suggests that, the GRB events 
exist in short and long duration classes, which are separated at 2 seconds [1]. In general, the 
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short duration bursts are harder and long duration bursts are softer. In the hardness ratio 
vs. duration plot, the two classes are seen to distribute in distinct domains [1, 2], where the 
hardness ratio HR 32 concerned is defined as the ratio of the total counts in the 100 — 300 keV 
and 50 — 100 keV energy range (later, the hardness ratio is defined as the ratio of the fluence 
in the 100 — 300 keV range to the fluence in the 50 — 100 keV range [3]). When sources of 
both classes are combined, the hardness ratio is obviously correlated to the duration, while 
for each of the two classes alone the two quantities are not correlated at all [4]. This in turn 
strongly suggests the existence of the two classes of GRBs. The study of GRB classification is 
very important since different classes might have different progenitors. It was proposed that 
the long duration bursts are caused by the massive star collapsars [5—7] and shot duration 
bursts are produced in the event of binary neutron star or neutron star-black hole mergers 
[8, 9]. 

After the successful launch of the Swift satellite [10], a large number of evidences favor 
the two progenitor proposal for GRBs. Long duration bursts were found to be originated 
from star-forming regions in galaxies [11], and in some of these, supernovae were detected to 
accompany the bursts [12, 13]. On the other hand, short duration bursts were detected in 
regions with lower star-formation rates, with no evidence of supernovae to accompany them 
[14-16]. 

Although GRBs are known to belong to the two distinct classes, the membership of the 
classes is hard to establish due to the overlap of the duration distributions, which makes the 
classification of the bursts an unsettled issue. Authors of Ref. [17] showed that, the bimodal 
distribution of GRBs can be well accounted for by two overlapped lognormal distributions, 
which suggests that each of the two GRB populations is likely to form a single peak modality 
and there would be a sufficient number of bursts that are mis-classified by simply applying the 
criterion T 90 = 2s. Authors of Ref. [18] used a multivariate analysis to discriminate between 
distinct classes of GRBs, and from the third BATSE catalog they found that instead of two, 
there are three classes of GRBs, which triggered a debate lasting until the present time. 

In distinguishing the X-ray flashes (XRFs), X-ray rich bursts (XRRs) and conventional 
gamma-ray bursts (CGRBs), authors of Ref. [19] used the ratio between the X-ray fluence 
Sx and the gamma-ray fluence S 7 as a softness parameter to divide them. Their results seem 
to be quite satisfactory. We wonder, if the softness parameter can play a role in separating 
the short-hard and long-soft populations of GRBs from the BATSE catalog. This is the 
motivation for our analysis given below. 

The BATSE catalog has the largest sample of GRBs available until today. We employ 
this sample in the following analysis since the total number of GRBs in the sample is large 
enough for a meaningful statistical analysis. 
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In Section 2, we investigate if the duration distribution depends on the softness param- 
eter. Based on this investigation, we shortly explore a possible classification of GRBs with 
an empirical curve in Section 3. The case of Swift is considered in Section 4. Discussion and 
conclusions are presented in Section 5. 



2. Duration distributions of GRBs with different softness parameters 

Motivated by Ref. [19] 's work, we wonder if the softness of BATSE bursts can play a 
role in separating the short-hard and long-soft populations of GRBs. However, most of the 
BATSE bursts do not have X-ray fluxes, so, we cannot define the softness of the bursts with 
Sx/S-y. Therefore, we turn to consider the ratio between the fluence in the lowest energy 
band and that in higher energy bands, as adopted in Ref. [20] to distinguish them. This 
parameter is no longer applicable to distinguish XRFs, XRRs and CGRBs. Here, we define 
the softness parameter of BATSE bursts by 

where f\, f'2, f'3, and are the fluences of the first (20 — 50 keV), second (50 — 100 keV), 
third (100 - 300 keV), and fourth (> 300 keV) BATSE channels, respectively. 

As hinted by Ref. [19, 20], we roughly divide the bursts into several groups according to 
the softness parameter SP. To be comparable to the conventional duration distribution plot, 
we divide the duration range from 0.01 s to 2000 s in the logarithm format into 30 bins. We 
require that the maximum count of each group should be no less than 30 to keep a statistical 
significance. For the softest group, we require that if there exist some short duration bursts 
in this group, the corresponding count should be noticeable. In this way, we get four groups 
of bursts. The distributions of durations (T 90 ) of these groups are displayed in Figure 1, 
where T 90 of a burst is the time interval during which 90% of the total observed counts have 
been detected. We find that the duration histogram of each group bears a bimodality. The 
positions of the dip lying between the two peaks of the bimodality of different groups differ 
significantly, where the softer the group, the smaller value of the duration position of the 
dip. If separating the bursts into short and long duration classes, then in different ranges 
of the softness parameter the classes occupy significantly different percentage of the total 
counts. The figure shows that the long duration bursts occupy a very large percentage of 
the softest population (SP > 0.1), and they cover a large range of durations, starting from 
~ Is to ~ 1000s. In contrast, short duration bursts occupy a very small percentage of the 
softest population, and they cover a very small range of durations. The situation changes 
for harder groups: the percentage of long duration bursts becomes smaller and they cover 
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Fig. I. — Duration distributions of the BATSE bursts with different softness parameters. 
Range of the adopted softness parameter is mentioned at the upper left corner for each panel. 
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smaller ranges of durations, while the percentage of short duration bursts becomes larger 
and they cover larger range of durations. 

In the conventional classification of GRBs, a dip around 2 s in the duration histogram 
plot in Figure la of Ref. [1] was proposed. For a more precise measurement, the authors fitted 
a quadratic function between the two peaks in the histogram and determined its minimum 
to be T90 = (1.2 ± 0.4)s. As an empirical analysis, we adopt this quadratic function fitting 
method to locate the dip position from the plots in Figure 1. In Table 1, dip positions and the 
numbers of long and short duration bursts divided by this feature for various softness groups 
are listed. In the first column, the softness parameter range of each group is presented. The 
mean of the softness parameter calculated with all the individual SP values of the bursts of 
the corresponding group is presented in the second column. Positions of the dip are given in 
the third column. The numbers of short and long duration bursts separated by the dip are 
presented in the fourth and fifth columns respectively. 



3. Empirical curve roughly setting apart two populations of BATSE bursts 

A relation between the softness parameter and the dip position of the bi-modal duration 
distribution of BATSE bursts hinted by Figure 1 can be evaluated from the data of Table 
1. In Figure 2, we display the result of a power-law analysis on the data. Fitting the data 
with a power-law function yields (by performing a Spearman correlation analysis with the 
ORIGIN software): 

SP = (0.100 ± 0.028)T 90 {0 - 85±0 - 18) . (2) 

The following are the statistical results of a linear correlation analysis between logSP and 
logT 90 : the correlation coefficient, -0.958; the number of data points, 4; the probability of 
rejecting the null hypotheses, 0.0417. Owing to the limit of the total counts, we have only 4 
data points in Figure 2. The statistical analysis performed here is thus not robust (we hence 
regard this fitting curve as an empirical function — see the statement below). Note that one 
can divide the BATSE catalog into more groups to have more data points. As a result, one 
might find difficulties to locate the dip suggested in Figure 1 due to data fluctuations. 

In Figure 3, we display the plot of the softness parameter vs. the duration of BATSE 
bursts, where the function of Equation (2) is also plotted. The figure shows that there exist 
two populations of bursts which are clustered in two distinct domains in the plot. Between 
the two clusters there is a "gap" with relatively sparse counts. The function of Equation (2), 
which is the fitting curve in Figure 2, is seen to pass through the "gap", roughly separating 
the two populations. As the figure shows, most sources of the population below the fitting 
curve are the conventional short bursts (the duration is smaller than 2 s), and most sources 
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Table 1. Parameters of different softness groups. 
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a The number of sort duration bursts with their duration being 
no less than the T90 value presented in the third column. 

b The number of long duration bursts with their duration being 
larger than the T 90 value presented in the third column. 
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Fig. 2. — Relation between the softness parameter and the position of the duration distri- 
bution dip drawn from the BATSE catalog. 
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Fig. 3. — Plot of the softness parameter vs. the duration of BATSE bursts, where the solid 
line is the function of Equation (2). 
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Fig. 4. — Duration distributions of the two populations divided by the curve of Equation 
(2), where the thick line represents the short-hard population while the thin line denotes the 
long-soft population drawn from the BATSE catalogue. 

of the population above the fitting curve are the conventional long bursts (the duration is 
larger than 2 s). 

We find that: a) the number of the population below the fitting curve is 511, and among 
them there are 433 bursts (84.7% of the population) with their duration satisfying T 90 < 2s; 
b) the number of the population above the fitting curve is 1461, and among them there 
are 1418 bursts (97.1% of the population) with their duration satisfying T90 > 2s. The 
duration distributions of the two populations divided by the fitting curve is shown in Figure 
4. Remind that the fitting curve represents the relation between the dip positions of the 
duration distribution of bursts and the softness parameter. According to this analysis (refer 
to Figures 1, 3 and 4), we regard the fitting curve as an empirical curve roughly setting apart 
the short-hard and long-soft populations of BATSE GRBs in the softness vs. duration plot. 
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4. In the case of Swift 

GRBs observed by Swift are becoming more and more important due to the excellent 
ability of Swift to detect the afterglows of the objects. The number of bursts collected by 
Swift has been steadily growing up since the launch of the satellite [10]. Would the statistical 
analysis performed above be applicable to the current Swift data set? 

To answer this question, we search in literature and collect the Swift data which contain 
the values of the duration and the four Swift channel fluences. The data are available in 
Ref. [21] where the first Swift BAT GRB catalog is presented. From Ref. [21] we get 222 
bursts with their T 90 as well as f l (15 - 25 keV), f 2 (25 - 50 keV), / 3 (50 - 100 keV), and 
fi (100 — 150 keV) being available. Compared with Table 1 we find that this number of 
bursts is too small to be used to perform the same analysis. This explains why we focus our 
attention on the BATSE catalog instead of the Swift catalog. 

Even though the number is so small, we manage to divide the bursts into two groups 
according to their softness parameters to check their duration distributions. The same 
definition of the softness parameter is adopted, although the energy channels of Swift are 
not the same as those of BATSE (also, T90 would not be exactly the same as that measured 
in BATSE, since the energy range of bursts would play a role in determining the quantity). 
Those Swift bursts with SP > 0.1 are included in the softer group, whilst the others belong 
to the harder group. Duration distributions of the two groups are shown in Figure 5. We 
observe that, although the statistical significance is much less than that in the BATSE 
catalog due to the small number of bursts, the same trend is observed in the figure: the 
duration histogram of each group is likely to bear a bimodality; the positions of the dip 
lying between the two peaks of the bimodality of the two groups are different, where the 
softer the group, the smaller value of the duration position of the dip; when separating the 
bursts of each group into short and long duration classes according to the apparent dip, 
then any of the two classes occupies significantly different percentage of the total count in 
different groups, e.g., the softer the group, the larger percentage of the long duration bursts. 
Comparing Figure 1 with Figure 5 one can find that the adopted statistical analysis is hard 
to be applied to the current Swift data due to the limited number of bursts. 

To check the trend in other ways, let us ignore the statistical significance and fit a 
quadratic function between the two peaks in the histogram of each panel of Figure 5 and 
determine the minimum of the curve as the position of the dip. The result is displayed in 
Figure 6, where as a comparison, the BATSE result shown in Figure 2 is also presented. The 
same trend of the empirical function is maintained if we rely on the two Swift data points 
in the figure. However, we would like to emphasize that one should not take the two data 
points so serious since the numbers adopted in Figure 5 are quite small so that the statistical 
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Fig. 5. — Duration distributions of Swift bursts with different softness parameters. The 
range of the adopted softness parameter is mentioned at the upper left corner of each panel. 



1 r 



Q- 



0.01 - 



0.1 r 




1E-3 



-1.0 



-0.5 



0.0 



0.5 



1.0 



1.5 



2.0 



log "U s ) 



Fig. 6. — Relation between the softness parameter and the position of the duration distri- 
bution dip drawn from the Swift catalog (open circles). Other symbols are the same as they 
are in Figgure 2. 

result is not so robust. In addition, we adopt the same definition of the softness parameter 
to analyze the BATSE and Swift catalogs but this definition would correspond to at least 
a slightly different quantities since the energy ranges of the corresponding channels are not 
the same. While this is not an unsolvable problem currently, but the number of bursts is. 
We hope that the growing Swift burst number would provide a reliable analysis in the near 
future. 



Motivated by Ref. [19] 's work, we investigate if the softness of bursts plays a role in the 
conventional classification of GRBs. We employ the BATSE catalog and define the softness 
parameter as the ratio of the fluence in the first channel to the sum of the fluences in the 
second, third and fourth channels. The duration distributions of different groups of BATSE 
GRBs associated with distinct softness parameters are explored. From the analysis we get 
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Discussion and conclusions 
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Fig. 7. — Distributions of the conventional hardness ratio (the thick line) and the softness 
parameter (the thin line) of the bursts in the BATSE catalog. 

an empirical curve that can roughly set apart short and long duration bursts. The analysis 
shows that the conventional classification of GRBs with the duration of bursts is influenced 
by the softness of the objects. 

As is known, the hardness ratio is an active factor to tell the spectral difference between 
the short and long duration bursts of GRBs. Why do we adopt the softness parameter 
instead of the conventional hardness ratio to investigate the issue? As mentioned above, it 
is Ref. [19] 's work on dividing XRFs, XRRs, and CGRBs with another softness parameter 
(S x /Sj) that motivates us for the exploration of the softness parameter as a possible factor 
active in the classification of GRBs. While the role of the conventional hardness ratio in 
GRB classification has been well studied for a long time, the role of the softness parameter 
has not been sufficiently explored. Another reason for doing so is illustrated in Figure 7, 
where distributions of both the conventional hardness ratio and the softness parameter are 
presented. We find that, while the conventional hardness ratio spans about one magnitude 
interval, the softness parameter covers approximately two magnitude range. This suggests 
that, if both the conventional hardness ratio and the softness parameter can be adopted as 
a parameter to distinguish different types of burst, then the latter must be more sensitive 
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than the former, since the difference of the latter is more easily noticeable. 

What would one get in the conventional hardness ratio vs. duration plot, if the bursts 
of the BATSE catalog are divided by the empirical curve, the function of Equation (2)? 
The result is shown in Figure 8 where we find that two distinct populations are clustered 
in the conventional domains of the short and long duration classes. The well-known corre- 
lation properties between the hardness ratio and duration of the two classical populations 
are maintained. That is, for any of the two populations, the hardness ratio is not at all 
correlated with the duration, while the two elements are obviously correlated when the two 
populations are combined. This is not surprising, since 84.7% bursts of the population below 
the fitting curve are the conventional short duration bursts (T90 < 2s), and 97.1% bursts of 
the population above the fitting curve are the conventional long duration bursts (T 90 > 2s). 
This together with Figure 4 show that dividing BATSE bursts with the empirical curve or 
with the duration criterion T 90 = 2s would not give rise to a statistical difference (or, the 
statistical properties of the two resulting populations would not be significantly different for 
the two criterions). This in turn suggests that the modification to the conventional method 
(e.g., by considering different softness groups instead of a single sample, or by adopting the 
empirical curve instead of the Tg = 2s criterion) would not lead to a noticeable physical 
discrepancy between the two populations. 

As Figure 8 shows, the two populations divided by the empirical curve heavily overlap 
in the hardness ratio vs. duration plot. This together with the overlapping in Figures 3 and 
4 imply that, the two well-known populations are unlikely to be sharply separated merely 
by quantities such as the duration, the hardness ratio and the softness parameter, or the 
relations between them. To distinguish the intrinsic short-hard and long-soft populations, 
other supplemental criterions are needed. Among them, the most desirable one is the location 
of bursts in the host galaxies, as observable in some Swift bursts. 

One might notice that the softness parameter adopted here is in fact an inverted BATSE 
gamma-ray hardness ratio HR 2 u/i = (/2+/3+/4) / fi- Therefore, the fact that the dip existed 
in the bimodality in the duration distribution of GRBs shifts with the softness parameter 
(see Figure 1) is equivalent to the fact that the dip position shifts with the hardness ratio 

HR234/1- 

What is the difference between our analysis and the conventional one? The conventional 
investigation considers a single sample from a catalog whilst our study considers different 
groups of softness from the same catalog. The conventional one yields a single dip position of 
the histogram of the duration distribution while our analysis produces several dip positions 
which are found to vary with the softness. The statistical methods adopted to distinguish 
the long and short duration classes in both cases are exactly the same, but the results 
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Fig. 8. — Hardness ratio vs. duration plot for the short duration bursts (filled circles in the 
upper panel; open circles in the lower panel) and long duration bursts (filled circles in the 
lower panel; open circles in the upper panel) divided by the function of Equation (2) from 
the BATSE catalog. 
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are different. Our analysis suggests that the conventional classification should be modified 
by setting apart the two populations in different softness groups, which would be more 
reasonable than doing this with a single sample. Another possible modification is to separate 
the two populations in the softness vs. duration plot with an empirical function as long as 
the function is reliable enough (when the total number of bursts is much larger than the 
current available one, the function so obtained would be more robust in terms of statistics, 
e.g., in that case one would get more data points rather than only 4 in Figure 2 under the 
same requirements). 

Assume that we take the empirical function as the criterion to classify GRBs. If so, is 
the new classification statistically better than the old one? As discussed above, as much as 
84.7% of the newly classified short-hard bursts are the conventional short duration bursts and 
as much as 97.1% of the newly classified long-soft bursts are the conventional long duration 
bursts. The new method is only a modification to the conventional one but not a totally 
different one. Therefore we cannot expect an entire difference between bursts identified by 
the two kinds of classification. However, if we believe that each of the two populations 
follows a lognormal distribution of the duration of bursts, as suggested in Ref. [17], then 
we can check the statistical difference of the two classifications in this light. Presented in 
Figure 9 is the fitting result of a lognormal distribution curve to the duration distributions of 
BATSE bursts for the two classifications. The statistical improvement from the new method 
is obviously observable. 

As Figure 9 shows, compared with the short duration class, the long duration class is 
more poorly fitted by a lognormal distribution. This indicates that if the empirical function 
is used as a criterion to separate the BATSE bursts, the resulting long-soft population might 
contain a "third" subclass. As mentioned above, about one decade ago, authors of Ref. [18] 
used a multivariate analysis to discriminate between distinct classes of GRBs and found that 
there exist three classes of GRBs instead of only two. The debate of the existence of the 
"third class" has then been triggered and lasting until the present time [22—28]. We suspect, 
if a detailed analysis based on what is shown in Figure 9 is helpful in ending the debate, and 
hence hope to see such an investigation in the near future (for example, to explore a fit to 
the duration distribution of the long-soft population with the method of the superposition 
of two lognormal distributions, as that adopted in Ref. [25]). 

What result would a different number of intervals of the softness parameter produce? 
Obviously, larger intervals would lead to a smaller number (smaller than 4) of data points 
in Figure 2, which will not be considered here. As mentioned previously, smaller intervals 
would give rise to more data points in the figure, but the statistical significance will be 
harmed. For example, when we divide the BATSE catalog into 7 instead of 4 groups, the 
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Fig. 9. — Duration distributions (the thin line) of BATSE bursts fitted by a lognormal 
distribution curve (the thick line) for the conventional short (the left bottom panel; the 
reduced chi-square for the fit is 85.6) and long (the right bottom panel; the reduced chi- 
square for the fit is 305) duration classes, and for the newly classified short-hard (the left 
top panel; the reduced chi-square for the fit is 35.8) and long-soft (the right top panel; the 
reduced chi-square for the fit is 227) populations by the empirical function of Equation (2), 
respectively. 
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Fig. 10. — Duration distributions of seven groups of the BATSE bursts divided according to 
the softness parameters. Range of the adopted softness parameter is mentioned at the upper 
left corner of each panel. 
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requirement that the maximum count of each group should be no less than 30 would no 
longer be satisfied (this is the reason why we prefer to divide the BATSE catalog into 4 
groups). As a result, more obvious fluctuation of data would be expected. The result is 
shown in Figure 10. Compared with Figure 1 we find that data fluctuation is indeed more 
obvious in Figure 10. Following the same method, we get an empirical curve from Figure 
10, which is SP = (0.104 ± 0.019)T 90 ( ' ' 95±0 ' 11 ' ) . Parameters of this curve are in agreement 
with Equation (2) (within the error bars, values of the fitting parameters obtained in both 
cases are the same), indicating that our analysis does not rely strongly on the choice of the 
number of intervals of the softness parameter. 

According to the above analysis, we reach the following conclusions, a) When divid- 
ing BATSE bursts into several groups according to their softness parameters, the duration 
distribution of each group bears a bimodality histogram, b) Between the two peaks of the 
bimodality exhibits a feature of dip that sets apart two subgroups, shorter duration sub- 
group and longer duration subgroup, c) The dip position shifts with the softness parameter, 
where the softer the group, the smaller value of the duration position of the dip. d) For 
groups with different softness parameters, the short duration subgroup occupy significantly 
different percentage of the total count, where the softer the group, the smaller percentage of 
the total count it occupies (for the long duration subgroup, the conclusion is the opposite), 
e) Deduced from the relation between the dip position and the softness parameter we find 
an empirical function that can roughly separate the short-hard and long-soft populations. 

As a primary investigation, we find a statistical difference between the duration distribu- 
tions of a single sample and multi-group samples from the BATSE catalog. We do not know 
what causes this difference. Nor are there any models that have ever predicted this. Lack- 
ing the knowledge of the causing, the physical interpretation associated with the proposed 
modification of the classification is currently unavailable. In the same way, the phenomenon 
of overlapping is currently not able to be explained. Other independent investigations are 
needed to understand the new findings as well as the overlapping phenomenon. 
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